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ABSTRACT 


Experiments, under different conditions, were conducted 
to evaluate the effect of condensate inundation on condensation 
heat transfer in tube bundles of marine condensers. Five 15.9 
mm. (5/8 in) nominal outside diameter, smooth stainless steel 
tubes were used in a vertical row to simulate an actual 
condenser. Tubes were located in an equilateral triangular 
array with a spacing-to-diameter ratio of 1.5. 

Heat transfer performance was determined for each tube in 
a bundle. Data was taken for condensing steam on the outside 
of each tube at about 21 KPa (3 psia) and at about 101 PKa 
(14.7 psia). Each tube was cooled by water on the inside at 
velocities of 0.78 to 6.22 m/sec (2.56 ft/sec to 20.42 ft/sec). 
The overall heat transfer coefficient was determined directly 
from experimental data. The inside and outside heat transfer 
coefficients were determined using the Wilson plot technique. 

Observation of condensate flow showed lateral droplet 
motion along the first three tubes in portions of the con- 
denser under all conditions tested. Side drainage occurred 
only over the third and fourth tubes at a condensation 
pressure of about 21 KPa. The dominate mode of the flow at 
101 KPa condensation pressure was gravitational flow. Outside 


heat transfer coefficients were higher than expected under all 





conditions when compared to the Nusselt theory. The reason 
for this is possibly due to secondary vapor flow. Recommendations 


to improve validations are provided. 
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te oN RODUC TION 


A. BACKGROUND INFORMATION 

Recent improvements in turbine machinery and boiler design 
have broughtabout an increase in horsepower to weight ratio of 
marine propulsion systems. However there has been no similar 
improvement, in practice, in the steam plant condenser size. 

Practical marine steam condenser design is based almost 
exclusively upon two documents. The Heat Exchange Institute 
(HET) standards for steam surface condensers [1] and the 
Standards of Tubular Exchange Manufacturers Association 
(TEMA) [2]. These standards have proven to be reliable. But 
it is evident that the resulting condensers are significantly 
overdesigned. Briefly, Search [3] has shown that heat 
transfer enhancement methods could decrease condenser space 
to weight ratio, thereby establishing new design criteria 
for marine condensers. 

It has been the objective of past research by Eshleman [4] 
to investigate the effect of condensate inundation on heat 
transfer in a horizontal tube bundle. But, his research 
concentrated on designing, constructing and validating the 
test facility. 

Since the publication of Nusselt's well-known theoretical 
paper [5] on film condensation many theoretical and 


experimental studies have occurred. 


LS 





The evaluation of the heat transfer coefficient from 
the vapor condensing on the outer surface of horizontal 
tubes is based mainly on Nusselt's theoretical formulas. 

Nusselt performed the derivation and found that the 
average heat transfer coefficient AL for n tubes located 


below one another is, 


3 4 
.. Pe (Pe-P,,) Degks g 
See) (a) 
nH~D. (TL -T) 
where Ke : Thermal conductivity of film (W/m- C 

Des ; Average tube wall temperature ( C) 

Pe Film density (kg/m?) 

P., Vapor density (kg/m?) 


Up : Dynamic viscosity of film (kg/m~sec) 
g : Acceleration of gravity (m/sec*) 


For the uppermost tube, equation (1) becomes 


3 4 





h ‘ = 0.725 
eee ee tee ) 

Hence, the relationship between the average heat transfer 

coefficient for a horizontal tube bundle (h.) One Loring Of 

n tubes and the mean heat transfer coefficient (Aya, ) for the 


uppermost tube may be found as, 


Rhee on = nA 
Nu (a) 
The above equations were derived by Nusselt assuming: 
1. Condensation of a saturated vapor at negligible 


velocity, 
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2. Laminar flow of the condensate film in a continuous 
sheet from one tube to the next at a constant 
temperature difference for all tubes in the bank. 

In actual condensers, vapor moves at a fairly high speed 
over a considerable part of its path. Under changing turbine 
speeds, steam velocity 1s not negligible. But attempts to 
evaluate analytically the effect of vapor velocity have not 
been successful. The vapor velocity causes friction between 
the vapor and the condensate film. With downward flow of 
vapor and condensate (as in our experiment) the frictional 
forces are added to the force of gravity. Consequently as 
the film velocity increases, the thickness of the film 
decreases and the coefficient of heat transfer from vapor 
to wall increases. 

Additonally, the condensate does not flow down from each 
tube in a continuous film. Rather it forms drops or streams. 
Condensate dripping on a tube from above splits around the 
tube but does not flow axially. The thickness of the film 
caused by the condensate coming down from above is thus 
confined to the place where the drops and streams descend, 
which causes local disturbances in film flow. Briefly, the 
true nature of condensate flow on and between tubes differs 
from the flows assumed by Nusselt. 

Recently, Eissenberg [6] noted that the condensate droplets 
formed on the tubes strike anywhere on the lower half of the 
tubes below. His experimental results gave heat transfer 


coefficients well above those predicted by Nusselt. 
i 





Eissenberg proposed a side drainage model, and 


formulated that, 


ee ncn eos 42 a (4) 


>t 


As a matter of fact, it can be said that in most cases 
experimental data has been fit to various modified forms of 


equation (5) below: 


h /h =n” (5) 
mmere 0.074sS0.223. 

Generally, the effects of condensate inundation and vapor 
velocity are described separately from each other. Actually, 
they occur simultaneously and their combined effect is 
complex. Fujii [7] correlated the effect of inundation and 
vapor shear by using experimental data of Nobbs and Mayhew 
[8, 9]. The data for in-line tube banks resulted in 
equations (6) and (7).. 


Nu? = 10.74 Reo’? ae 
m L 


Equation (6) above is an experimental equation for a tube 


without inundation where 


° 


Nu: Nusselt number, (hod. /ky,) for pure steam, also 


without inundation. 
Re, : Two phase Reynolds number, (Udo /v,) 
iaeemakl Gonduetivity o£ liquid 
dad. : Tube outside diameter 
Jeeeekinemakie Viscosity of Liquid 


and U_ : Vapor velocity 
16 


Nu_/Nuy is computed for the data with inundated tubes 


and equation (7) is derived. 


Nu_/Nu- = ( Re, /2x10° 0.072 (w/w) °° Gary) 
Nu >: Mean Nusselt number for a tube 
w : Rate of inundation falling onto a tube 
w- : Rate of condensation of a tube correspond- - 


; O 
SLIGeMES Tu, 


For the staggered tube bank, a correlation was not 
obtained by Fujii.Fujii determined that for the staggered 
tube bank, the inundation effect was smaller than that for 
the in-line bank.This would follow from examination of 
Eissenberg's side drainage model. 

Chisholm [10] combined the developments of Berman and 
Tumanov [11] , and Fuks [12] into one eauation to evaluate the 
effects of condensate inundation, vapor velocity and non- 
condensable gas in tube bundles. Chisholm has ociven the 


following formula for heat flux, 


4 
ag ae a _ 
where 
-0.07 
n-1 
ee 
K> Lg re i 
eno 725 | ——— eee) (100095 Re °7 YN) (94 
h 2 m 
u W 
Ore Ce eon 1 


na 


See 





Equation (9) is for the downward flow of the vapor. 


T. : Temperature of condensate surface (°C) 

T., : Temperature of outer tube wall (°C) 

K : Thermal conductivity of condensate (W/m-°C) 
L : Specific latent heat (J/kg) 

ve: Absolute viscosity of condensate (kg/m-sec) 


V_ : Specific volume of condensate (m?/kg) 
: Rate of condensation of rth tube row (kg/sec) 


Re. : Reynolds number of vapor-gas mixture. 


In the above equation, 


o D 
stat o 
— Kk. Coron) 
where ae Heat transfer coefficient across the condensate 


LMM state le Vapor Condltion (Ww/m*- Cc) 


Experimental data resulting from varying inundation 
rates and vapor velocities is sufficiently scattered to 
Suggest that no existing correlation fits all the available 
data (some of the available data is related to enhanced 
tubes). This is because there are so many variables that 
affect marine condenser performance, some of which have to 
be discovered (e.g., non-condensable gas effects, pitch-to- 
diameter ratio of condenser tubes, direction of vapor flow, 
and the effects due to ship motion in three dimension). 
Research is presently being conducted in the United Kingdom 
aimed at producing a model of condensation suitable for 
sophisticated condenser performance calculations. However, 


the results have yet to be published. 
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eB OBJECTIVES OF THIS WORK 

In order to evaluate the effect of jnundation on a bank 
of enhanced tubes, it was necessary to establish data ona 
bank of smooth tubes as a standard of comparison. The 
objectives of this work were therefore: (1) to establish 
experimental data on a bank of smooth tubes as a standard, 
and (2) to compare the established data to theoretical 


predictions. 
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IIT. EXPERIMENTAL APPARATUS 


foe LNTRODUCTION 

The existing test facility was designed by Beck [11] and 
built and tested by Pence [12]. Major modifications were made 
to the original apparatus by Eshleman [4]. During this work, 
some minor modifications were made and these changes are given 


meerables IIT and III for each run. 


eee OLkEAM SYSTEM 

The steam system is shown in Figure 2. The supply of steam 
is locally generated and supplied to the building which houses 
the experimental apparatus. The steam is provided by means of 
a 19.05 mm. diameter line and a steam inlet valve (MS-2). A 
compound pressure gage is located just prior to the steam 
separator which monitors the supply pressure as it is adjusted 
by (MS-2). Following the steam separator, a line strainer 
provides additional protection from contamination. After the 
strainer, the steam proceeds through a 31.75 mm. diameter line 
which provides for two possible steam paths. The primary path 
for system operation is via the throttling valve (MS-3), 
through a desuperheater and into the test condenser. Inside 
the condenser, the steam is condensed on the test tubes. The 
steam which is not condensed proceeds via the vapor outlet on 


the test condenser to the secondary condenser. The secondary 
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steam flow path is used to accomplish system stabilization 

mmeing Startup and to control the mass flow rate of steam to 

the condenser during operation. Steam proceeds via (MS-4) 
directly to the secondary condenser. All steam lines except 

the primary path downstream of (MS-3) were insulated with 25.4 mn. 


thick fiberglass insulation. 


fee LEST CONDENSER 

The test condenser is shown from various views in Figures 3 
and 4. Steam enters via the top and proceeds over the baffle 
separators and through a flow straightener, which is covered 
with three layers of 150 mesh screen, to the tube bundle. The 
condensate collects at the bottom of the condenser and flows 
out one of the two 12.7 mm. diameter openings at either end of 
the condenser to the hotwell where it can be collected and 
measured. 

Three separate viewing windows each 203 mm. by 140 mm. by 
17.7 mm. and made of pyrex plate glass had been installed to 
provide maximum viewing of the active tubes. 

The tube sheet arrangement is shown in Figure 4. The tubes 
were arranged in a typical naval condenser spacing-to-diameter 
Metro Of 1.5. They were 15.9 mm. OD, 1.14 mm. thick, 304 
stainless steel tubes that had cooling water passing through 
them. Although typical naval condenser tubes are made of 90-10 
copper-nickel, the choice of 304 stainless steel was based on 


"on hand" stock and the fact that the principles of inundation 
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do not depend on the tube material although perhaps the heat 
mlax may change due to different wall resistance. The remaining 
half tubes were made of 15.9 mm. OD aluminum bar stock and were 
fastened by screws to the outside wall of the steam flow path. 
This arrangement was selected to best simulate the steam flow 
conditions in a section of an actual condenser. The five test 
tubes are sSingularly removable. The top tube can be replaced 
by a 304 stainless steel porous tube which could simulate various 
condensate inundation rates. 

The test condenser was insulated with a 25.4 mm. thick sheet 


er Armorflex insulation. 


mee CONDENSATE SYSTEM 

The condensate system 1S shown in Figure 5. The condenser 
hotwell collects the condensate from the test tubes, while the 
secondary condenser hotwell collects the condensate from the 
secondary condenser. Valve (C-l1) allows the isolation of the 
test condenser hotwell for condensate measurement. Since house 
steam was used as the steam supply system, the condensate 
collected in the hotwells is pumped back to the house system 
by the condensate pump via valve (C-3). The condensate lines 
were insulated with 19.1 mm. thick Johns-Manville Aerotube 


maisulation. 
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E. COOLING WATER SYSTEM 

The cooling water system is shown in Figure 6. The water 
used was normal house water which had been passed through a 
water softener on the way to the supply tank. The water is 
eemeed from the supply tank by a 5 HP electric driven pump. 
Tt is routed to the flowmeter header via 51 mm. OD plastic 
pipe. The flow of cooling water for each test tube is then 
individually controlled and measured by it's own rotometer. 
Each rotometer allowed a maximum flow rate of 70.4 LPM. The 
heated cooling water, after passing through the test section, 
Was piped back to a supply tank. A separate system pumped this 
heated water through a filter and cooling tower returning the 
cooled water to the supply tank in an effort to maintain a 
constant cooling water inlet temperature. 

After leaving the rotometers, the system piping was reduced 
to 15.9 mm., ensuring a distance of 1m. ahead of the test 
section to ensure a hydrodynamically fully developed velocity 


profile while passing through the test section. 


fee soeCONDARY SYSTEMS 
1. Vacuum System 
The vacuum in the test condenser and secondary condenser 
was maintained by a mechanical vacuum pump and a vacuum regulator 
which induces air into the system. The vacuum pump takes a 
suction from the secondary condenser hotwell which is connected 


to the test condenser hotwell via discharge piping. A cold 


ZS 





trap at the inlet of the vacuum pump forces incoming vapor 
to pass over a system of refrigerated copper coils. This 
removes steam and entrained water from the vacuum line prevent- 
ing moisture contamination of the vacuum pump oil. The vacuum 
pump outlet is vented through a roof exhaust fan to avoid a 
health hazard from breathing any oil vapor exhausted by the 
pump. A schematic diagram of this system can be found in 
Figure 7. 
2. Desuperheater System 

The desuperheater removes sensible heat from the super- 
heated steam by injecting water at about 25°C via the existing 
feedwater system through valve (DS-1) and a rotometer. The 
desuperheater is a 267 mm. diameter stainless steel can, 457 mn. 
high, having four nozzles inserted equidistant around the 
circumference of the inner top of the can. The nozzles are a 
fan type and are positioned such that the spray is downward at 
a 45° angle to allow for better mixing. A collection tank is 
located on the bottom of the desuperheater to allow for drainage 


of condensate. This system can be isolated by valve (DS-2) 


G. INSTRUMENTATION 
1. Flow rate 
a. Steam velocity was determined by calculation 
m Vv 


Ss Cli) 
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where 


Th = Mass flow rate of condensate (kg/sec) 
= Whe, 
A. = Cross sectional flow area a 
and y = Specific volume of vapor are ep 


b. Cooling water flow rate was measured individually 
for the five separate tubes. Each flow rate was determined 
by a rotometer with a capacity of 70.4 LPM (18.6 GPM). The 
calibration procedure used was identical to that listed in 
Appendix A of Ref. [14]. 
2. Pressure 
Two different pressure sensing devices were used during 
experimentation. They were a Bourdon tube pressure gauge which 
measured steam pressure and an absolute pressure transducer 
coupled with a 760 mm. mercury manometer which was used to 
measure test condenser pressure. 
3. Temperature 
Stainless steel sheathed , copper-constantan thermo- 
couples were used as the primary temperature monitoring devices. 
Figure 3 shows the location of five vapor thermocouples. The 
remaining 30 themocouples of this type were located as shown in 
Figure 6, Six on each tube, two measuring cooling water inlet 
temperature and four measuring water outlet temperature. 
Calibration procedures of the thermocouples were identical to 


those listed in Appendix A of Ref. [16]. 
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4, Data Collection and Display 
An Autodata collection system was utilized to record 
and display the temperatures in degrees celsius obtained from 
the primary stainless steel thermocouples. Table I lists the 


channel numbers and location of these devices. 
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ie Aree lMEN ca PROCEDURES 


A. OPERATING PROCEDURES 
i «= PueGparation of condenser tubes 

Prior to any run, each tube was properly prepared to 
ensure filmwise condensation. The cleaning procedure for 
stainless steel tubes is listed in Appendix A. 

ee system operation and steady state conditions 

The basic operating instructions developed by Pence 
[14] and modified by Reilly [15] were used. The only differ- 
ence in the procedure as listed in Appendix B of Ref. 16 was 
that instead of one cooling water flowmeter to adjust, the 
experimenter had five to set as desired. 

In general it took about three hours from initial 
light-off until steady-state conditions were established. 
The parameters used in determining steady-state conditions 
were cooling water inlet temperature and steam inlet 
temperature. If the cooling water inlet did not vary more 
than +0.6 °C/hr and the steam temperature did not vary more 
than +0.3 °C/min, steady-state was considered achieved. 

The time for the system to stabilize between changes 
in cooling water flow rate during the Wilson plot technique 
was approximately twenty minutes. This time increment is 
suspect as other investigators waited about one hour for 


stabilization between changes, especially for atmospheric runs. 
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It must be pointed out however that the amount of time required 
to collect data over five tubes in a system that can't be shut 
down and repeated the next day prohibits the greater time 
increment between data points for the Wilson plot. 

The general set up for the data taken in this 
research was a steam velocity of between approximately 1.2 
m/sec and 1.4 m/sec, steam temperature of 62 °C or 75 °C for 
21KPa condensation conditions and 100 °C for atmospheric 


(101 KPa) runs. 


B. DATA REDUCTION PROCEDURE 

The raw data collected for each tube for run 10 can be 
found in the Tables beginning on page 38. 

Appendix B, the sample calculations, is a complete listing 
of the equations used to evaluate the data. Appendix C is 
a derivation of the probable error in the data reduction 
equations, followed by a sample error analysis for tube 
number 1 at 40 percent flow, run 10. 

The following standard heat transfer equations were used 
to reduce the raw data into a form that can be used for 
evaluation. 

1. Overall heat transfer coefficient (U,,) 

The method employed to arrive at the overall heat 
transfer coefficient is straightforward and similar to that 
employed by many researchers in the past. The heat transfer 


rate to the cooling water is given by 
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ce ex = T .) ie) 


The heat transfer rate from the steam is given by, 


= 1 Caer + = 
2 Moon Ico. Vv ae ee Econ cat nde | a 


The heat transfer rate can also be found from the overall 


heat transfer coefficient by 


Q = U_ A, (LMTD) ( 14) 


where 


ie eS Ge 
V Gi V eo 


LMTD = 





( 15 ) 
Lie Tee Tt — Toe) 
V ei V ce 


After combining equations (12), (14), and (15) it can be 


Beund that 
m a ty _ err 
es ( 16 ) 
Mm Vv CoO 


fp orbected overall heat transfer coefficient (U,) 





CE) 





where So is the wall resistance corresponding to different 


tube materials and may be given by 
A in(x /x,) 


so @ fies) 


2. nek 
W ae 


where Ly. is the length of the tube. 


Equation (17) allows for the comparison of tubes of different 
Materials for the same steam and cooling water condition 


within the test condenser. 
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3. Inside heat transfer coefficient (hj) 


need, 
ey 
Nu = = 0.036 Pe. pri/3 (ee) ee ( 19 ) 


k 





Equation (19) was selected because both the Dittus-Boelter 
and Sieder-Tate relationships which are commonly used 
assume a fully developed velocity, as well as a fully devel- 
oped thermal profile. In this research, it was suspected 
that, although the velocity profile was believed to be fully 
developed, the thermal profile was not fully developed. 
When an L/D ratio of 57.6 is used in equaiton (19) a constant 
of 0.029 results. This was validated by computing the aver- 
age of all the tube constants obtained as a result of the 
Wilson plot technique. Wilson plots for each tube for run 10 
can be found in Figures 8 through 12. 
PoLelde heat transter coefficient (h,) 

The outside heat transfer coefficient is the parameter 

that is used to compare results of each tube in the bundle and 


is given by : 


i (2205) 
(1/uU_) - RR - (D/D;h, ) 
Two very important assumptions were made in using this 
equation. 
a. The resistance due to fouling was equal to Zero. 
This assumption is supported by the fact that the tubes were 


new, chemically cleaned and smooth. Also, treated soft water 


was used as the cooling medium. 
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b. The resistance due to non-condensable gases was 
equal to zero. This assumption was supported by the fact 
that the system was tested for air-tightness and found to be 
secure. In addition, it was believed that the velocity of 
steam passing through the test section was sufficiently large 
to keep the system purged of any non-condensables that might 


collect in the test section. 


C. DATA REDUCTION COMPUTER PROGRAM 

Reilly [15] developed the existing program in Fortran 
Language. His program had been translated into Basic Langu- 
age for use with the HP 9845 computer by Eshleman [4] during 
his-work. Ultimately, this program with minor modifications, 
can be used in an integrated system between the Autodata Nine 
data collector and the HP 9845 computer. This will allow 
automatic data input with real time data output for the 
experimenter. During this work, the existing computer 
program of Eshleman [4] for reduction of data was used with 


the HP 9845 computer. 
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ihe ercoUbto AND DEISCUSS.1ON 

The experiments were done by using two different conden- 
Sation pressures to establish experimental data on a bank of 
smooth tubes.Condensation pressure for runs 1 through 5 and 
9 through 12 was maintained at about 3 psia. For runs 6, 7, 
and 8 the pressure was maintained at atmospheric conditions. 
Runs 8, 9, and 2 were repetitions of runs 7, 5, and l 
respectively.Experimental conditions are given in Tables II 
emo Lit for each run. 

In equation (16), the T term stands for vavor temperature. 
Large differences in the heat transfer coefficients were 
obtained depending on whether the actual vapor temperature or 
the saturation temperature was used in equation (16). It is 
worthwhile to note that the outside heat transfer coefficients 
using the actual, superheated steam temperatures were 52 per- 
cent lower than those calculated using the saturation 
temperature.The percent change of oe values was 79 percent. 

To overcome this interesting result, it was decided that the 
T, term in equation (16) must be the saturation temperature 
instead of actual vapor temperature. For comparison U. and 
a values are tabulated in Table XIV for runs 10 (saturated) 
and 11 (superheat) in which T., was actual vapor temperature. 
Saturation temperature was therefore used in equation (16) 


for all runs except those reported in Table XIV. 
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The runs, conducted at 3 psia condensation pressure, gave 
unexpected results. Outside heat transfer coefficients for 
the five test tubes showed significant fluctuations. The 
values of h, decreased for the first three tubes and then 
increased. Generally, for all runs conducted at 3 psia, the 
heat transfer coefficients followed this same pattern. The 
cause for this behavior may be due to several phenomena. 

During observation of the condensate flow pattern on the 
first three tubes, there was evidence of lateral droplet 
migration. This migration was due presumably to a deflection 
of these tubes or due to axial flow of vapor. This causes 
a non-uniform heat transfer rate across the length of the 
tube because of a decrease in film thickness at some locations 
and then an increase in film thickness in other locations. 
This thickening of the film on the tubes would result ina 
lower average heat transfer coefficient than expected. 
Secondly, steam appeared to concentrate at the bottom of the 
test condenser, and may have caused cross flow around the 
lower tubes. This may have been caused by a system resist- 
ance to the flow of excess steam which is not condensed in 
the test condenser. Circulating flow of the excess steam 
which could not easily leave the test condenser may have 
caused side drainage of the condensate. It was observed 
that the condensate flow path on the fourth and fifth tubes 
was toward the observation window. Because of these two 


reasons, the outside heat transfer coefficients for the last 


a3 





two tubes would be high. As an example, Table XXII lists 
Values for run 10, and these are plotted in Figure 13 along 
with two theoretical results. 

Figures 8 through 12 are the Wilson plots that assist in 
determining the constant in equation (19). As an example, 
run 10 was chosen and the results are plotted in these 
Figures. All tubes yield good linear plots with slopes which 
Beeovide constants of 0.028, 0.032, 0.032, 0.028 and 0.030 
respectively. The expected linear plots were obtained for 
all the other runs as well as run 10. The data reduction 
program gives the option of using the constant solved for 
via the Wilson plot technique or inputing one of the user's 
own choosing. In this work, for all runs, 0.029 was used as 
input for all the tubes to determine the inside heat transfer 
coefficient which in turn was then used in the determination 
of the outside heat transfer coefficient. Outside heat 
transfer coefficients for all runs are tabulated in Tables XV 
meeeough XVIII. 

The ratio of h/hy as listed in Tables XIX through XXII 
for all runs, was determined by taking the average outside 
heat transfer coefficient Ay for n tubes and averaging them, 
then dividing by the outside heat transfer coefficient of 
the first tube (hj). The results of run 7, as listed in 
Table XXI, are plotted on Figure 14 along with the theoretical 
equations of Nusselt and Eissenberg. Based on the observa- 


tions of condensate flow at atmospheric pressure, the data 
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for all tubes was expected to fall closer to the Nusselt 
curve due to the presence of gravity dominated flow, but 
certainly not below it. 

The experimental study on the effect of the vapor velocity 
upon condensation was performed both at 3 psia and 14.7 psia 
condensation pressure. Runs 9 and 10 at 3 psia, and runs 6 
and 7 at atmospheric pressure were conducted at different 
vapor velocity to determine the effect of vapor velocity. 
Vapor velocities for runs 10 and 7 were slightly higher than 
runs 9 and 6 respectively. The experimental result for the 
3 psia pressure case wasS unexpected. This result may be due 
mo, the BErect of saturated steam. At atmospheric pressure, 
for the higher vapor velocity, the heat transfer coefficient 
was higher as expected. 

From examination of the results of runs 7, 8 and 10, 1l 
1t may be said that at higher condensation pressure, the heat 
transfer coefficient is higher. 

Each repeated run was within the uncertainty range of the 
Original. The conclusion is that the data obtained for ten 
runs is sufficient to make comparison between results. 

Side pieces and baffles were installed inside the test 
condenser to prevent axial flow of the vapor. However droplet 
Migration still resulted in both cases. From the comparison 
of the run 1l and 12, it may be concluded that side pieces 
generally did not have any effect on the heat transfer 


Goetficient. Heat transfer coefficients of runs ll and 12 
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are within the uncertainty range of each other. The differ- 
ence between heat transfer coefficients of runs 4 and 10 is 
due to effect of baffles. Outside heat transfer coefficients 
are higher at run 4 because of the baffles that form a path 


for the condensate. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The experimental data obtained lead to the following 


conclusions: 


1 


There is evidence of secondary flow of steam within 

the test condenser which is suspected to have influenced 
the results. 

Equation (16) should not be used with superheated steam 
temperatures. 

The measured ratio of h/hy is lower than expected at 


atmospheric runs. 


The following recommendations are provided: 


Ihe 


Improve the steam flow path to ensure a uniform downward 
profile through the tube bundle. This can be accomplished 
by re-design of the steam inlet section. 

Prevent auxiliary system resistance against the flow 

of uncondensed steam when working at vacuum. This can 

be done by either conducting experiments at atmospheric 
pressure or by re-deSigning the auxilary system piping 
with large diameter. 

Instrument the top tube with thermocouples to measure 

i in order to check the outside head transfer coefficient 
against the Nusselt theory. 

Measure the vapor pressure between tubes to get more 


accurate calculations. 


By 
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TABLE I 
LOCATION OF STAINLESS STEEL SHEATHED COPPER-CONSTANTAN THERMOCOUPLES 
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eee 0 iT 
RUN CONDITIONS FOR RUN 1 THROUGH 7 


Condensation pressure 3 psila. ,Superheat 


a : 
Seca. Low vapor velocity. 
9 Samemes cum No. i 
Condensation pressure 3 psSia.,Superheat 
3 


steam,Low vapor velocity,Steam gage 


has been changed. 


5 Same as) hUneNOe; 4+ 
Condensation pressure atmospheric, 
Doe eeeamoeeaml, Low Vapor Velocity, 


Baffles have been removed. 


Condensation pressure atmospheric, 


Z 
Saturated steam,High vapor velocity. 





Condensation pressure 3 psia.,Saturated 

4 steam,Low vapor velocity,Baffles and 
Vacuum control valve have been installed, 
Throttle valve has been changed. 
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TABLE III 
RUN CONDITIONS FOR RUN 8 THROUGH 12 


_ Same as Run No. 7 


Condensation pressure 3 psia,Saturated 




















steam,Low vapor velocity. 


Condensation pressure 3 psila. ,Saturated 


Stcamewten vapor velocity. 


Condensation pressure 3 psia.,Superheat 


steam,High vapor velocity. 


Condensation pressure 3 psia. ,Superheat 
steam,High vapor velocity,Side pieces 


have been installed. 
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TABLE ALV 


CALCULATED Un AND h., VALUES FOR RUN 10 AND 11 BY USING ACTUAL 
VAPOR TEMPERATURE 
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OUTSIDE HEAT TRANSFER COEFFICIENTS FOR RUN 1, 2 AND 3 
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APPENDIX A 


TUBE CLEANING PROCEDURE 


To ensure filmwise condensation, the condenser tubes had 


to be prepared. Surfaces of the tubes were cleaned to insure 


proper wetting characteristics and to insure that all deposits 


were removed. Stainless steel tubes were prepared in accordance 


with the procedure given in Newton [17]. The steps in this 


cleaning procedure are as follows: 


i. 
Dri 


Prepare an Alconox detergent solution and heat to 90 °C. 
Apply this solution to the surface of the tubes. 

Drain and rinse the tubes with distilled water. 

Spray with alcohol. 

Rinse with distilled water. 

Spray with acetone 


Rinse with distilled water. 
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APPENDIX B 


SAMPLE CALCULATIONS 


The following is an example of how the data reduction 


program calculates the results. Tube number one at 40 percent 


flow rate of cooling water of run 10 was selected for this 


analysis. This same tube and flow rate was used for the error 


analysis in Appendix C. 


Input parameters 

Tube Outside Diameter (D.) 

Tube Inside Diameter (D; ) 

Tube Length ae 

Outside Nominal Surface Area (A) 
Wall Resistance (R) 

Cooling Water Inlet Temperature ae 
Cooling Water Outlet Temperature lee) 
Average Cooling Water Temperature (Thor TEE) 


Steam saturation temperature a) 


Gallons Per Minute of Cooling Water (GPM) 


Section 1. Water Properties 


yu (MHUW) 


O05 3875 mM. 
02014) im. 
0.9144 m. 


ienercmric 


> mMnek/H 


5.72x10 
2 eG 
B27 ae 
Bone ns0s.4 | K 
73.375 °C 


7.68 GPM 


eae EXp | ['2-a0829258) (303.2)« (2644. 2184) /(303.4)] 


emo 25206 } 


u Se ole x Me) kg/m-sec 


va 





kK (KW) = 0.5565919 + (0.002174417) (30.2) - (0.70127 x 10”) (30.2)° 


= (280914 x 10~°)(30.2)> 


K = 0.615858 W/m-°C 


P(RHO) = 1004.44434 - UG2075365) (30). 2) 


~ ( 0.0023913147) (30.2) 


sos. 436 kg/m 


Co mee) = 4¢-237/955 — (0.0018553514 ) (30.2) 


+ (1.3948314 x 10>) (30.2) 


Cp = 4-195 kI/ kg-°C 


ih (MFRCW) = LPM x RHO x 1.67 Sales 


where LPM=GPM x 3.78533 
a = (29.0713344) (998.436) (1.67x10 >) 
nh = 0-484732 kg/sec. 


Prandtl Number (Pr) 


uc / Meee 10) % 4.1945x 10°) / (0.615858) 


er 


Pr = 35-3315 


Section 2. Data Reduction 


2 
imecoOoling water velocity ce) no ee 
= (4 x 0.484732) / [ (998. 436) Tt (0.0141) “| 
C =) 3.1092 m/sec 
cw 


UA 





2. Mass flow rate per unit area (G) 


‘ Z 


G = (998.436) (3.1092) 


G = 3104.3372 kg/m--sec 


3. Reynolds Number (Re) 


Re =D,G / u = (0.0141 x 3104.3373)/ (7.828 x 10 *) 


Re = 55916.14 


4. Overall Heat transfer Coefficient (U,,) 


m C Tee 1 ne 
U = p in Vv ez 
1 A ie Su sk 
V 





n co 
= (0.484732) (4.1945 x 10°) Pie SoS 27 ad 
0.0456 ee Se eo. 


EuaOMace W/m = °C 


9. Corrected Overall Heat Transfer Coefficient (U .) 


Cc il 1 ~5 
7 eee Miesmigge =» > - * CC® 


TaBORGCnwe) m= °C 


U 
Cc 


6. Wilson Plot Parameters (X,Y) 


(a) Ordinate 


ae 
U ~ 5169.408 
Te 
ese x 10” ee Gl / w 
(b) Abscissa 
ee eee Ck ~ 
Pee ee Cay (ses ane) 7 


x = 9-114 x 10> 


Ve 





Determination of constant 


D 
Oo 


Mk 





oa 


Where M = Slope returned by linear regression subroutine 


M = 0.9265 
_ 0.015875 

(0.9265) (0.615858) 
era 0028 


Inside Heat Transfer Coefficient (he) 








h.D. 
f= ——— = 0.036 oe pa = (ly D wee 
k Oo 
where 
i 0.9144 _ 
Oo  Omolss75 — ae 
a = s (0.029) ere oA 
ys De 
a 
_ 0.615858 0.8 1/3 
ee 029) (5916.14) (5.3315) 
h, = 13898.64 HUeeo 
Outside Heat Transfer Coefficient (h) 
= 2 
oO 
a 
a eer 
U W Dee 
n aye ae 
h = s 
1 we 0.015875 
5169-400 >°/* * 10° ~~ ~(OTpiai) (13898.64) 
h_ = 18103.22 W/ rue 
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APPENDIX C 


ERROR ANALYSIS 


The basic equations used in this section are reproduced 
from Reilly [15]. The general form of the Kline and 
McClintock [18] "second order" equation is used to compute 
the probable error in the results. For some resultant, R, 
which is a function of primary variables Xy1 Kor ceeeee , Xr 


the probable error in R, 6R is given by: 


2 Zz 2 
sR = SR sy in OR sy eee ea 
6X 2 6X 2 6X n 
1 2 n 


where 5X1, SX5, ecu 6x, 1s the probable error in each of 
the measured variables. 


1. Uncertainty in overall heat transfer coefficient, U. 
6U, 6A, Z 6c, 2 6m \ 2 
ia 2) es oe ee oe : 

n n p m 


ewe ys 
V CL 


(T -T .)(T -T _) Ln ——— 
ae we 2. ee Ce ne ee 
2 2) 
L OTK | ; ST oo 
os - ~ = 
Coe.) Ln Ty Toi (Po—r..) Ln Ty Fed 
V CL 7 -T V CoO 7 -T 
Vewce Ve CO 


The following values are assigned to the variables 


iv 





éA_ = £0.0001 m@ 


6C,, = 40.0042 Kd/kg-°C 


6m = £0.01 kg/sec 


= 9 
= ° 
oT aa asin sxe 
— 0 


Ber tube No:]1 at 40 percent of 10. run: 


su, MOGI On0042 \~ 0.01 \? 
i 0.0456 4.1945 Oe4Sa 732 


2 
KOR) i= >) 
(567s) 4407675) tn(l.1229) 


1/2 
2 


fe eel all 0.1 
(40.675) In(l.1225) |; (45.675) Ln (1.1229) 


O03 7 


2. Uncertainty in inside heat transfer coefficient, h, 
The probable error in the inside heat transfer 


coefficient is given by: 
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a sx \° + ae _ {0.8 6Re : i cee se cbr 
h; k D; Re Este 


where 
ok = +0.001 W/m-°C 
6D; = +0.001m 
ae — +0210 


6C = +0.001 


P. 2 2 
i 0.001 0.001 2 
a ( ats ) ( 2a) + (0.8x0.073) 
2 2 
0.10 0.001 

7 (0.333 Ate) + ($2598 ) 

Sonaos 
h,, Tube No:1_40%_10.Run=13899 + 1362 W/m*-°C 


oy 


Y 





3. Uncertainty in outside heat transfer coefficient, ae 
The probable error in the outside heat transfer 


coefficient is given by: 

















2 
ey et D. 
; Ya SO, 7 Bw 7 Rr? a, Sw 7 DFR; 
eel: 
Sh pe 
De aoe ey 
Dh, oe 
ad 
l D 
fa Se O 
Un w «DA; 
where su, = 0.037 
Un 
aR = 2.86x10~° m-=°C/W 
sh 
ae = 0.098 
ee 
ms 25000 
Si = 5.524x10 7° m“-°C/W 
ee dnd 
I =e 7 
Sh, 0a 2.86x10 
_ (5169. 408) (5.524107) es s24= lor 
i 
2 2 
, | (810x107?) (0.098) 


(Ses2d<lOn>) 
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Sno 
h 
Oo: 


= 0.200 


h., Tube No:1-408-10.Run = 1810343621 W/m*-°C 
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